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Abstract 

We discuss phases of gauge theories in the holographic context, and formulate a criterion 
for the existence of a Higgs phase, where the gauge redundancy is "spontaneously broken", 
in purely bulk language. This condition, the existence of a finite tension solitonic string 
representing a narrow magnetic flux tube, is necessary for a bulk theory to be interpreted as 
a Higgs phase of a boundary gauge theory. We demonstrate the existence of such solitons in 
both top-down and bottom-up examples of holographic theories. In particular, we numerically 
construct new solitonic solutions in AdS black hole background, which are used to demonstrate 
that the bulk theory models a superconductor, rather than a superfluid. The criterion we find 
is expected to be useful in finding holographic duals of color superconducting phases of gauge 
theories at finite density. 
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1 Introduction and Conclusions 



Holographic methods are becoming a standard tool in analyzing quantum field theories and 
revealing physics which would be strongly coupled using other, more conventional descriptions. 
Most relevant to this note, this basic tool has been utilized in contexts in which we expect 
"spontaneous breaking of gauge invariance," such as holographic superconductivity (for a 
review see [l]) or color superconducting phases in QCD (for recent attempts to model such 
phases see (2|[3]). 

Of course, the expression "gauge symmetry" and its breaking is a misnomer, or more 
precisely relies on specific classical limit for its definition. In a specific weak coupling limit 
it makes sense to speak of gauge redundancies as approximate global symmetries and use the 
machinery and language of global symmetry breaking in this context. However, in an inherently 
non-perturbative context such as holographic dualities one needs to stick to more precise and 
gauge-invariant definitions. Such characterization of massive phases of gauge theories was 
given by 'tHooft [4||5], and we review this classification in section 2. 

This classification of gauge theory phases is gauge invariant and non-perturbative, relying 
on the response of the gauge theory vacuum to massive external sources. This could be best 
used in the holographic context whenever we have an idea of the gauge theoretic microscopic 
definition of the system, and use the holographic context merely to perform calculations in 
the strongly coupled regime. This situation is demonstrated in section 3, using one particu- 
larly simple such "top-down" context, namely that of the Coulomb branch of the maximally 
supersymmetric SU{N) theory in four dimensions. We demonstrate that the phase structure 
of the theory is manifested in certain geometrical features of the bulk theory which reproduce 
the expected results. 

The purpose of this exercise is to extract a purely bulk criterion for the existence of a Higgs 
phase interpretation of the theory, which we can then use in situations where the microscopic 
definition of the bulk theory is less well-understood. Indeed, we see that the expected behavior 
of the 'tHooft loop operator in the Higgs phase implies the existence of certain type of solitonic 
strings localized in the IR region of the bulk theory, representing a narrow magnetic flux tube 
in the boundary theor}|^ In the holographic context, this can be taken as the definition of such 
phases, since it implies much of the phenomenology we associate with the Higgs mechanism. 

Since our criterion depends only on the bulk geometry, it is ideal in the bottom-up approach 
to holographic duality, where the microscopic definition of the theory is lacking. In section 
4 we use this criterion to demonstrate that the Marolf-Ross prescription [6^ (see also j7|) for 
obtaining boundary dynamical gauge fields, can be used in the broken phase of the theory. We 
show that in that broken phase, those "dynamical" boundary conditions lead to holographic 
superconductivity, rather than superfiuidity. 

We do so by constructing new solitonic solutions in AdS^ black hole background (in the 
probe limit), by solving numerically the bulk equation of motion - a set of coupled non-linear 
partial differential equation. Section 4 is devoted to setting up the equations and boundary 

'*The precise interpretation of this flux tube depends on the microscopic interpretation of the theory, and in 
particular on the UV region of the geometry. 
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conditions, and describing the properties of the solutions. Essential to our solutions is the use 
of dynamical boundary conditions for the gauge fields (introduced in |6] ) , which as we will show 
are necessary for obtaining finite energy solutions, corresponding to superconducting vortices. 
We describe the dependence of their energy on the temperature, below the critical temperature 
for the onset of superconductivity. As expected, the solitons we construct disappear above the 
critical temperature. 

We are hopeful that the criterion discussed here, and the role it plays in models of holo- 
graphic superconductivity, will assist in formulating the problem of holographic color super- 
conductivity, and in constructing holographic models along the lines of ^ . We hope to return 
to this problem, one of the original motivations of the present note, in the near future. 



2 Characterization of Gauge Theory Phases 

In [4|[5] 'tHooft introduced a classification of phases of gauge theory based on its response 
to electric and magnetic sources. For the characterization to be a precise definition of the 
associated phases, we restrict ourselves for now to theories with gauge group SU{N) /Zj\f, 
such as gauge theories based on unitary groups in which all matter fields are in the adjoint 
representation. In such theories the centre of the gauge group Z^r is a global symmetry which 
aids in providing order and disorder parameters to characterize the different phases. 
The response of the theory to electric sources is measured by the Wilson loop 

W{C) = Tr{e'fc^) (2.1) 

where we take the trace in the fundamental representation. The curve C is taken to represent 
the worldline of two static external sources separated by distance L , and the Wilson line then 
computes the static potential between these sources. 

The response to magnetic sources is similarly represented by a 'tHooft loop T{C), which 
plays a role of a disorder parameter in the theory. The 'tHooft loop operator is defined in 
the path integral language as an integral over all gauge field configurations with a prescribed 
singularity along the curve C. The singularity represents the presence of an external magnetic 
sources. For the curve C which represents the worldline of two well-separated static sources, 
this operator probes the theory in a way which is similar to the Wilson loop. Indeed, as 
is well-known, these two observables are exchanged under electric-magnetic duality (see for 
example |8] or section 10 of Witten's lectures in [9|). 

We can then distinguish the different phaseqjof gauge theories by following asymptotic 
behavior for large loops C: 

• Confinement W{C) ~ e-^(<^) and T{C) ~ e'^C^) 

• Higgs Phase W{C) ~ e^-^C^) and T{C) ~ e'^^*^) 

""This is not a complete classification of such phases. For example, there could be critical points and oblique 
confinement phases distinguished by the behavior of dyonic loop operators. We will not discuss such phases here. 
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We denote the area enclosed within the curve C by A{C) and the corresponding behavior 
of the loop operator is called the area law. This encodes the linear potential between the 
corresponding (electric or magnetic) sources. The linear potential has an intuitive picture in 
terms of the existence of flux tubes connecting the sources (confining strings) which in turn 
exist because the corresponding (electric or magnetic) flux lines emanating from the sources 
form narrow flux tubes and do not spread (the Messiner effect). Similarly, the length of the 
curve C is denoted by L{C), and the corresponding behavior for the loop operator is called the 
perimeter law. Such behavior encodes the fact that the flelds generated by the corresponding 
source are short ranged (screened) and influence only the close vicinity of the source location. 

In the Coulomb phase, or in a conformal fleld theory, the behavior of both the Wilson and 
'tHooft operators is dictated by conformal invariance. For the loops corresponding to static 
sources separated by distance L, we have the behavior 

W{C) ~ T{C) ~ e-^ = e-^^(^) (2.2) 

where T is a large time cutoff, and a is a constant (which can depend on coupling constants 
of the theory). Note that while formally this is classifled as a perimeter law for both Wilson 
and 'tHooft loops, the behavior of the static potential V{L) in a massive (screened) phase is 

L 

different V{L) ~ e where Lq is the screening length. 

The prescription of calculating the Wilson and 'tHooft loops in AdS/CFT is simple and 
well-knowi][^ The expectation value of Wilson and 'tHooft loop respectively, in the funda- 
mental representation and when working in the saddle point approximation, is of the form 
e~^ . The action S is the minimal action of the worldsheet of fundamental string or D-strings 
respectively, in a configuration which end on the prescribed curve C on the boundary. As quan- 
tum operators the Wilson and 'tHooft loops obey an interesting algebra which constraints the 
possible phases of gauge theory, which was discussed in the context of AdS / CFT by Witten 
(section 5 of [lO]). 

The qualitative behavior of the Wilson loops in confining theories is also well-known. The 
electric fiux tube connecting external sources is mapped into a string worldsheet dipping into 
the bulk. In the confining phase the electric flux lines are confined to narrow flux tubes. The 
dual statement is that the string worldsheet localizes in the bulk radial direction, oftentimes 
for clear geometrical reasons (e.g. the "end" of the IR geometry in some sense). In the next 
two sections we provide an analogous statement for magnetic flux tubes in holographic theories 
in the Higgs phas^ 
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We use in our holographic discussion the BPS loops, the so-called Wilson-Maldacena loop and its magnetic 



dual 11 12 . The asymptotic behavior for the loops we consider is unaffected by the presence of the scalar fields. 



For suggestions on calculating the Wilson loop itself see 



^For previous discussion of these flux tubes, see 14 . 
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3 Top Down Model 



Consider k flat probe D3 branes in AdS^ x located at r = v (in Poincare coordinates), 
and smeared over the sphere S^. Here v is proportional to the VEV of the adjoint Higgs field 
giving rise to the Higgs mechanism in the = 4 SYM theory. This corresponds to the pattern 
of symmetry breaking SU{N) — )• SU{N — k) x SU{k), in the large limit, while k is kept 
finit^ In this example we have the power of large as an organizing principle, and we'll 
see that it aids us in separating the effects of symmetry breaking on the electric and magnetic 
loop operators. 



Electric Flux Lines 

We are mainly interested in magnetic flux tubes, but we start with a brief discussion of the 
Wilson loop. In the broken phase, with the above breaking pattern, the static potential 
between electric sources in the fundamental representation is schematically of the form 



a , k e 



—cvL 



V{L) = - + b 



L N L 

where a, b, c are constants. The leading order potential is still Coloumb-like, but since k gauge 
bosons are now massive we have a correction involving exchange of those massive gauge 



bosons. This can be seen, for example, if we repeat the calculations of 15, 16 in the broken 
phase. 

In the bulk this modification can be explained simply, as follows. The Wilson loop calcu- 
lation corresponds, in the saddle point approximation, to finding the area of a fundamental 
string worldsheet whose boundary ends on the prescribed curve C. The leading order term 
in the expansion contributing to (3.3) corresponds to the calculation in pure AdS 
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The form of the leading correction in (3.3) suggests a modification of the action of the same 
saddle point. 

The required modification arises when considering the worldvolume theory on the probe 
D3 branes. Consider the worldsheet of the fundamental string for well separated electric 
sources, in the broken phase, a situation which is depicted in figure (|3]). The leading order 
contribution for the Wilson line corresponds to the worldsheet area, and corrections come 
from interactions between the part of the worldvolume intersecting the probe branes (two lines 
on the probe branes, represented by two points in figure ([s])) . Since the radial fluctuations of 
the probe branes are massive - those correspond to the longitudinal modes of the W-bosons 
- it is easy to see that exchange of the massive scalar fields corresponding to these brane 
fluctuations reproduce the form of the leading correction in (|3 



*^This is conventionally called to Coulomb phase, and indeed the leading order interaction between electric sources 
will be Coulomb-like. Nevertheless we'll use the term Higgs or broken phase. 
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Probe Branes Boundary 



Figure 1: Wilson loop stretched between widely separated sources on the boundary. The leading 
order correction in the expansion comes from exchange of massive scalar representing the radial 
fluctuations of the probe branes. 

Magnetic Flux Lines 

In contrast to the calculation of the Wilson loop outlined above, the 'tHooft loop expectation 
value changes character from perimeter to area law, already in the leading order in the jj^ 
expansion. This corresponds to the existence of a new type of saddle point, rather than a 
modification of the action of the existing worldsheet. 

The new saddle point is similar to that of the Wilson loop in confining theories. Indeed, 
in such case the geometry of the bulk provides an IR cutoff, such as a soft or hard wall, or 
cap to the geometry. The area law is realized geometrically as the Wilson line for widely 
separated electric sources receives contributions predominantly from the vicinity of the IR 
geometry. This is the holographic dual to the statement that the flux lines connecting two 
electric sources do not spread out in the confining vacuum. 

1 In our case the magnetic dual to that statement cannot be explained in terms of the bulk 
geometry alone. Indeed, the area law for the Dl brane has to arise from differences between the 
worldvolume theory of such brane and that of a fundamental string (for example the different 
dilaton coupling [l4|). In our simple model this is easy to identify: in the presence of the probe 
branes the worldvolume of the Dl branes can take a detour through the probe D-branes which. 



Probe Branes Boundary 



Figure 2: World volume of Dl brane stretched between widely separated magnetic sources on the 
boundary. The area law for the 'tHooft loop results from the existence of a string-like object 
localized in the radial direction. In this model such object can be represented as soliton on the 
worldvolume of the probe branes, drawn in a thick red line along the worldvolume of the probe 
brane. 

for widely separated magnetic sources, will minimize the action. This is due to the fact that on 
the worldvolume of the probe branes, the Dl brane can be transformed into a solitonic string of 
finite tension. Therefore asymptotically in such separation, the minimum action configuration 
would be the one depicted in Q in which the D-string worldvolume stretches mostly along 
the worldvolume of the probe branes. Note that this is qualitatively similar to the Wilson line 
in a confining theory, in that the radial location of the loop is stabilized at some fixed radial 
location for widely separated sources. 

The interpretation of the flux lines in this simple example depends in various ways on 
understanding the full gauge-gravity duality. In particular, we have used large scaling 
to distinguish electric from magnetic flux tubes, and correspondingly confinement from the 
Higgs mechanism. Furthermore, the microscopic interpretation of the theory helped identify 
the type of charges available in the gauge theory, and which can be connected by those flux 
tubes. Nevertheless, we have identified a necessary condition for the existence of Higgs phase 
interpretation of the theory: the bulk spacetime should support a finite tension solitonic object 
which is approximately localized in the radial direction. The existence of this object, dual to 
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a narrow magnetic flux tube, is necessary for the 'tHooft loop of the boundary theory to obey 
an area law. In the next section we demonstrate the existence of such solutions in a simple 
bottom-up model of holographic superconductivity. 



4 Application to Holographic Superconductivity 

In this section discuss a specific 2+1 dimensional bottom- up model of holographic supercon- 
ductivity [It]. As discussed above, an area law for the 'tHooft loop is guaranteed by a finite 
energy vortex solution of the bulk fields localized at a particular value of the radial direction, 
representing magnetic flux tube in the boundary theory. When we do not have a microscopic 
definition of the theory, we take the existence of such soliton as the definition of the Higgs 
phase in the bottom- up holographic context. We demonstrate below the existence of such 
finite energy solitons in the present context. 

Crucial to the analysis is the prescription given in [6] (see also [^) for obtaining dynamical 
gauge fields in the boundary theory, by requiring the bulk gauge fields to obey specific type of 
boundary conditions at the UV, which well refer to as "dynamical" boundary conditions. We 
show that with these boundary conditions the required vortex solutions exist, and furthermore 
have finite energy per unit length. This indicates that the model, in the broken phase, describes 
a genuine superconductor. 



The Model 

We work in the context of the bottom- up model of \17\. The action is: 



S 



1 

2k 



d^Xy 

6 

~~ L2 



R 



4 



I {df, - iqA^) 



where tit? < 0, and q is the charge of the scalar field. We work in the probe limit 
as: 
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q (X e 



A 



tx e 







(4.3) 
defined 

(4.4) 

In this limit the matter energy momentum tensor scales as and drops out of the Einstein 
equation, and the metric is unaffected by the matter fields, while the Maxwell and scalar 
equations remain unchanged. The background solution features an AdS Schwartschild black 
hole geometry which, for certain values of the thermodynamic variables (the chemical potential 
H, or equivalently the temperature T) develops a profile for the scalar condensate and the 
temporal component of the gauge field. This signals the onset of symmetry breaking below 
the critical temperature. 

We choose to work in cylindrical coordinates (t, p, 9, w) with the conformal boundary lo- 
cated at w = 0. Our background metric is then: 

r2 

-^(-f(w)dt'^ + /(wy^dw"^ + + p'^de^) (4.5) 



ds^ 
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In addition we make use of the scaling symmetries of our action to scale the horizon to = 1 
and we take the AdS radius of curvature to be L = 1. This fixes the metric function to be 
f{w) = \ — w. By dimensional analysis we expect all physical quantities to be proportional 
to the ratio of T / fx. In what follows we fix = 1 and examine the behavior of the various 
quantities as a function of T. 



Ansatz and boundary conditions 

We now discuss matter excitations to the homogeneous background. Guided by the known 
vortex solutions of the Abelian Higgs model in flat spacetime (reviewed in appendix A), we 
propose the following ansatz for the solutions we seek: 

^ {Aq{w, p),Q, Ae{w, p),Q) 
if) — )• il^{w, p) exp{in6) (4.6) 

where s is the topological number associated with the vortex solutiorj^ 

The equations of motion consist of the two Maxwell and one scalar equation: 



q^ApR^ ^ d,Ao , d^Ao 
2T~ + + + ^"'^o = 

«(^-(^.4-^).(4.^)c..iJ.%^ + ^.aJfl^O (4.7, 

where we rescaled the scalar field as ip{w,p) — )• wR{w,p), for reasons of numerical stability. 
It can be seen that in the probe limit described above there is a scaling symmetry of the 



equations (4.7), implying that if a solution is found for a given value of q it is known for all q 
via an appropriate rescaling of the fields. This property is convenient for numerical purposes 
as it allows us to choose a scale for the matter fields which is numerically tractable. 

We wish to solve our system of PDEs on the domain defined hy wq < w < and 
< /3 < oo, where wq is a UV cutoff. For the problem to be well posed we must choose self 
consistent boundary conditions which are also compatible with the bulk equations of motion. 
We choose the following boundary conditions on the four different segments of the boundary: 

• p — )• oo: In flat space it is known that the vortex fields decay exponentially towards 
asymptotic values for the gauge and scalar fields as p goes to infinity. Anticipating 
similar behavior in our case suggests that we impose a Neumann boundary conditions 
at some finite and large value, pcut, since in that region the solution should tend to the 



homogeneous ground state 
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^Of course, this number is not conserved in the full geometry, and indeed as well see that it "unwinds" as function 
of the bulk radial coordinate w. 

^'^Pcut is chosen such that our solutions vary by less than 0.01% if it is increased. 
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p — )• 0: To determine the boundary conditions at the vortex core we require that ah 
components of the bulk magnetic field be finite. The radial and transverse components of 
the magnetic field are given by = dpAg/p, and Bp = d^Ag/p, respectively. Finiteness 
of the radial component implies that dpAg — t- as p — t- 0. Regularity of the transverse 
component then restricts the Ag component to obey (in this limit) dwAg — )• 0. Therefore 
we conclude that Ag{w,p = 0) must be a constant. If we were to impose a Dirichlet 
conditions at the conformal boundary this would fix this constant to be zero. In our case 
we have a residual gauge freedonf^ consistent with the boundary conditions, which we 
use to set Agi^w, p = 0) = at the core of the soliton. 

w ^ wq : On the conformal boundary we impose Dirichlet conditions on the scalar 
and ^0 fields — the scalar field must be normalizable and niust asymptote to the 
chemical potential p. Crucially, we impose a specific Neumann boundary condition on 
the Ag field, demanding that its radial derivative vanishes at the conformal boundary. 
This corresponds to having a theory in which the boundary value of the bulk gauge field 



corresponds to a gauge field [6] in the boundary theor} We refer to these boundary 
conditions as the "dynamical" boundary conditions in what follows. 

w — >■ w+: Regulairty conditions at the horizon are necessary since the equations degen- 
erate there. Choosing the solutions which are regular at the horizon means that the 
coefficients of the divergent terms in a power series expansion of the equations near the 
horizon have to vanish. This prescription yields the following constraints in our case: 

Ao = 

q^Ag^ , '^qsAg _2 \ _ p , dpR , 



R + - - ^ - 3 - 3d^R +^+d'R = (4.8) 

\ P P P J P 

R^q (s - qAg) - 3d^R - ^ + d^Ag = 

P ^ 

We numerically solve the equations with these boundary conditions using successive over- 
relaxation (SOR) in the domain wq <w < Wj^ and < p < pcut^ where the truncation radius 
Pcut is large but finite. In this approach the equations are discretized on a lattice covering 
the domain of integration. We use a second order finite differencing approximation in which 
the derivatives are replaced with their finite differencing counterparts. An initial guess for the 
value of the scalar and gauge fields is then assigned to each grid point. Dirichlet boundary 
conditions are implemented by insisting that the initial values assigned to the fields at the 
boundary grid points are maintained throughout the relaxation procedure, whereas Neumann 
or Robin boundary conditions must be imposed after each iteration. This is done by using 
the discrete form of the derivative operators to update the boundary grid points based on the 



^"'^The dynamical boundary conditions at the conformal boundary allow for gauge transformations whose parameter 
is independent of the radial coordinate w. 

^^We choose to make dynamical only the component Aq of the gauge field, for the sake of simplicity, we do not 
expect the features of the solution to change much if is made dynamical as well, since it already has nearly 
vanishing radial derivative near the conformal f boundary, in all the solutions we are interested in. 
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values calculated for the interior points. The SOR algorithm then provides an iterative method 
of finding numerical solutions to this finite difference system to within a prescribed tolerance. 
Once the solutions are available other quantities of interest such as the energy density are 
calculated via insertion of these solutions into the suitably discretized action. Further details 
of our implementation are found in Appendix B. 



Free energy 

Before presenting the numerical solution and discussing its properties, we explain the reason 
we expect the energy (per unit length) to be finite in our case. The discussion parallels that 
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The Lagrangian density of the bulk fields is: 

L = aVVi?! _ ^2 (^(-Me)! + + r^) _ lvo^pf{d^Rf - w^pfRd^R - \w^p{d,Rf 

+^ + - ^^^^f^ - ^ (4.9) 

Since the resulting action diverges near the boundary, we regularize it by subtracting the 
action of the translationally invariant hairy black hole solutions from the on-shell vortex action. 
Such subtraction automatically removes the divergences which occur due to integration in the 
w direction. Therefore divergences, if they exist, can occur only as a result of p integration. In 
the region of large p the scalar and fields asymptote to their values in the translationally 
invariant ground state. Therefore the only terms in the (asymptotic) Lagrangian density to 
survive the subtraction procedure are: 



1 rpcut 
2Tr I dt I dw / dp 



RHs-qAef , w^fidy^Ae)^ 



+ 



2p 2p 



(4.10) 



where all fields are understood to be functions of w only. Here we have introduced the cutoff 
Pcut in order to regulate potential divergences in the p integration. 
If we now use the Aq equation of motion to make the substitution: 

+ (4.11) 

2 ^ ^ ^ 2q 2q ^ ^ 

and integrate by parts, using the fact that / vanishes on the horizon, we obtain the logarith- 
mically divergent term: 



Trlog ( — ] I dt( -w^ f Aedu,Ae \w=q + {-] I dwdw{up' fdwAe 







= vrlog (^^^ j dt {-w^fAed^Ae + -^fd^^A^ U=o (4.12) 

In integrating by parts we have introduced the length scale 7 which is a measure of the size of 
the vortex core. 
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This reasoning led the authors of [18] to conclude that their vortex solution is logarithmi- 
cally divergent, as expected from vortices in a superfluid. We see that if we instead consider 
dynamical boundary conditions for the Ag field, 9«,^e|to=o = 0, and provided an appropriate 
vortex solution exists, the coefficient of the logarithm will vanish. We see below that indeed 
such vortex solutions (whose profile significantly differs from the superfluid vortices found 
in [Is]) do exist and we calculate their finite energy (per unit length). This demonstrates that 
our model describes a genuine superconductor. 



R(w,p) Ag(w,p) 




p 



Figure 3: The matter field profiles at a temperature of ^ 0.89Tc In order to aid in visualization the 
background translationally invariant solution has been subtracted from the Aq gauge field. Note 
the asymptotic approach of the scalar and Aq fields to their translationally invariant profiles and 
the fact that Ag field is independent of the radial coordinate w, and asymptotes to - as p — oo. 
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The Solutions 



We are now ready to present our numerical solutions for the bulk fields and discuss their 
properties, including the profile of the energy density, both in the bulk and boundary theory. 
We also plot the temperature dependence of the free energy for a range of temperatures below 
the critical temperature. We leave discussion of our numerical solution to appendix B. 

The system has a critical temperature Tc, below which it is in the condensed phase (i.e. 
the scalar field develops a normalizable background). Below that critical temperature vortex 
solutions start appearing, in figure (3) we show the profile of the fields for a typical vortex 
solution. The form of the solutions may be understood as follows: far from the vortex core 
the fields tend to their homogeneous profiles and the PDEs reduce to ODEs. When solving 
these ODEs numerically one finds that the solution for the Ag field is a constant, given by 
|. Together with our previous discussion of the p — )• boundary conditions, this means that 
Ag asymptotes to a constant, independent of the radial coordinate, both as /) — t- and as 
p — )• oo. Since we are also demanding vanishing radial derivative at the conformal boundary, 
a reasonable guess is that the global Ag solution depends only on p, i.e. Ag{w,p) — )• Ag{p). 
This is indeed what we find numerically. As seen in (4.12) above, the asymptotic form of Ag 
is directly responsible for the finiteness of the vortex energy. 




Figure 4: Energy density profiles for the translationally invariant background and vortex solutions 
along a radial slice at p = 0.1, at a temperature of ~ 0.89Tc 

Once we obtain the numerical solutions for the matter fields, their (discretized) energy- 
momentum tensor is evaluated. In figure (4) we illustrate the radial profile of the energy 
densities of both the translationally invariant and vortex solutions. The subtraction of these 
two profiles gives us the regulated energy density, Tqq^, of the vortex solutions. We note that 
the bulk energy density of the vortex solution dips below that of the homogeneous ground 
state near the conformal boundary. Nevertheless, the boundary energy density of the vortex 
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(relative to the background) is everywhere positive. 

Since our solution is time independent, and we are working in the probe Umit, we can find 
the energy density of the boundary theory by simply preforming the radial integration of the 
bulk energy density with appropriate measure factor: 

eip) = -[ dr^jZ' - - E 0(^o7).. (4.13) 

i * 

where in the last expression we replaced the continuous coordinates p, w by their discrete 
counterparts Wj, used by our numerics. The boundary energy density is displayed in figure 
(5) for several values of the temperature. At low temperatures (relative to the critical temper- 
ature) one sees that the vortex energy profiles are sharply peaked near p = and that, as one 
approaches the critical temperature, they flatten and broaden as the vortices begin to disperse. 
The vortex solutions merge with the homogeneous background at the critical temperature Tc- 



j,,o* Boundary free energy density 




Figure 5: Free energy density for several values of the temperature. Notice the changes in the 
vortex profile as a function of temperature — at low temperatures it is peaked near the vortex 
core while as the temperature increases it tends to become wider and more diffuse, tending to the 
homogeneous background at the critical temperature T^. 

Figure (6) shows the decrease in the free energy as we approached the critical temperature 
(found from examining the translationally invariant condensate). It was however found to be 
difflcult to determine the precise behavior of the solution in the vicinity of the critical point. In 
this limit the energy profile becomes infinitely wide and shallow while the background energy 
density also tends to zero. The combination of these effects make the numerics intractable. 
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Free energy of vortex solutions near critical point 
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Figure 6: The free energy as a function of temperature below the critical temperature. We note 
that for sufficiently low temperatures the free energy increases in an approximately linear fashion. 

Conclusion 

In conclusion, in this section we constructed vortex solutions in the context of the holographic 
models of fl7\. These vortices signify the onset of local symmetry breaking. The imposition 
of the dynamical boundary conditions corresponds, via the prescription of |6], to a dual field 
theory with dynamical gauge field. This is evidenced, for example, by the fact that any 
boundary gauge transformation which is only a function of the boundary coordinates respects 
the dynamical boundary conditions on the bulk gauge field. We find that in the spontaneously 
broken phase, the symmetry breaking is manifested by the existence of bulk vortex solutions 
with the expected properties of superconducting vortices: there is no operator corresponding to 
a superfiuid current on the boundary, and the vortex boundary energy is finite. In contrast, as 
found at |18|, the imposition of Dirichlet boundary conditions leads to a theory which exhibits 
a global symmetry breaking and vortices with diverging energy, as expected in a superfiuid. 

We expect that the criteria developed here for characterizing local and global symmetry 
breaking will have applications in other bottom-up holographic models. In particular it would 
be interesting to explore the applicability of these techniques to models of finite density QCD 
and color superconductivity. 
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Appendix A: Vortex Solutions in Flat Spacetime 

The vortex solutions for the flat space Abelian Higgs model are well known. Here we provide 



a brief discussion following |19j, for further information see also 20 
The action of the Abelian Higgs model is: 

L = -\f^'^F^, + ]^{V,-^nv^-^) - ^A(|V^|2 - f2)2 

The constant F is proportional to the VEV of the charged scalar field ^ breaking the U{1) 
gauge symmetry. Finite energy configurations of the fields require that they obey the asymp- 
totic conditions: 

l^-I^Vo, V = {df,-^ - iqA^-^) ^ as x ^ oo 

Requiring cylindrical symmetry, the form of the scalar and gauge fields at infinity are con- 
strained to be: 

^'(r, 9) ■00 exp(ia((9)) 

i Snip 1 da 

^^-^ 1=—^' ^ r-^oo 

q ip qr ad 

The winding of the phase, a, at infinity is an integer, n, which is related to the quantised 
magnetic fiux through the plane orthogonal to the magnetic vortices: 

9 Z"^" e / e 

n = — / AardO = — i) A.dl = — x magnetic flux 

In order to see explicitly the localized nature of these vortices one is required to examine the 
equations of motion. Using the ansatz 74^ = ^0 = 0, Aq = A{r) and ^{r,0) = ip{r) exp{in6) 
we obtain the following equations: 



%lj{r) 



n 



qAir)]^ + Xiipirf - F"^' 
r ' 



r 
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The falloff conditions imply that ip{r) — )• tpQ and A[r) — )• — as r — )• oo. We may use this 

qr 

information to linearize the Maxwell equation for large r by setting -(/'(r) = ipQ. Solving the 
resulting equation yields: 

A{r) > 1 — 1= exp i-qFr) 

r-5>oo qr yjr 

As the scalar and gauge field approach their asymptotic configurations at large r, we expect 
the action to be dominated by the potential term. Therefore in order to find the asymptotic 
behavior of the scalar field we examine perturbations of the potential. Requiring that the first 
derivative of the potential to vanish fixes the minimum at V'o = f-- The fluctuations around 
this are of the form F'^\p(r) where p{r) is the deviation of the scalar field from ^o- Using this 

approximation for the potential in the scalar equation and setting Air) = — we obtain: 

qr 

il){r) )■ ipo + C2 exp {-VXFr) 

The localized nature of the vortex is evident from the exponential decay of the fields to their 
asymptotic values for large r. The full solution can be obtained by solving the equations 
numerically. 



Appendix B: Details of the Numerics 

We solve the equations with the boundary conditions listed in section |4] numerically using suc- 
cessive overralaxation (SOR) algorithm. To this end we discretize the equation on a lattice of 
finite mesh-size h covering the domain of integration, such that continuous spatial coordinates 
{p,w) are represented by discrete pairs (wi^pj), where 1 < i < Nw,l < j < Np are integers. 
We use a second order finite differencing approximation (FDA), where the derivatives are re- 
placed with their finite differencing counterparts, e.g. d^R — )• — Ri^ij) /2h, dpR — )• 
(-Rjj+i — i2jj-|_i)/2/i etc. Following discretization, we thus obtain finite difference equations, 
at every mesh point, for each field. We iteratively solve the entire system of algebraic equations 
using pointwise SOR starting with an initial guess for the fields, until a desired precision is 
achieved. Typically we initialize our scalar and ^0 gauge field with the values of the homoge- 
neous solution (found by solving the ODEs using shooting) and set the Ag field to its expected 
asymptotic value of n/q. Along the horizon an initial guess is made for the scalar and Ag 
fields which interpolates exponentially between the zero boundary condition at p = and the 
expected asymptotic values at the /o — )• 00 boundary. We use similar SOR parameters for all 
fields. These are calculated at each step via Chebyshev iteration. In this iteration the spectral 
radius of the Jacobi iteration is chosen, for simplicity, to be that of the Laplace equation with 



Dirichlet boundary conditions, see Section 19.5 of 21 for further details and the algorithm. 
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While Dirichlet boundary conditions are implemented by assigning the fields their initial 
values throughout the relaxation procedure, Neumann or Robin boundary conditions are up- 
dated after each iteration. We do this by using the backwards FDA derivative operators to 
update the boundary grid points based on the values calculated at the interior points. It turns 
out that at the horizon such a straightforward implementation of the regularity conditions 
(4.8) is numerically unstable. As these conditions relate the radial and tangential derivatives 
of the fields along the horizon they yield, upon discretization, a pair of coupled polynomial 
equations which relate the values of the fields at grid points in the near horizon region. At- 
tempting to solve these polynomial equations to update the values of the boundary grid points 
after each iteration resulted in instabilities, which we attribute to the fact that the linearized 
scalar equation near the horizon is ill-posed (the effective mass terms and the elliptic operator 
have the same sign). The physical reason of this instability can be traced to the fact that the 
effective scalar mass in the near horizon region violates the Breitenlohner-Freedman bound, so 
that it triggers an instability and formation of a condensate. 

We found that a stable implementation of the constraint equations (4.8) is to evaluate all 
terms in the equation, except for the radial derivative, on the line of grid points just before the 
horizon, and to then use the FDA form of the radial derivative to extrapolate to the values of 
the fields on the horizon. This approach is consistent with the bulk equations of motion and 
identical to implementing the desired constraint equations when the continuum limit is taken 
(i.e. the limit in which the step size is taken to zero). 

While our numerical lattice extends all the way from the horizon w = 1 to the conformal 
boundary w = 0, it covers only finite domain in the transverse direction < p < pcut- The 
truncation radius pcut is chosen such that our numerical solutions are altered by less than 
0.01% when pcut is increased. Typically we use pcut ~ 70 — 80. In addition, we checked that 
asymptotically our PDE solutions of the vortex configuration converged to the ODE solutions 
of the translationally invariant configuration at the transverse boundary to accuracies of 0.01% 
or higher. 

Finally we discuss convergence of our finite-differencing numerical solutions. The rate of 
convergence is assessed based on the assumption that in the continuum limit, when the grid-size 
tends to zero, the discrete solution on the mesh h, designated u^, approaches the continuum 
solution, u*, namely Uh = u^-\-0{h"'). The power n measures the rate of convergence. It can be 
calculated by running simulations with similar parameter settings on a sequence of meshes with 
decreasing mesh-spacings h,h/2 and /i/4, and computing n = log2{uh — Uh/2) / Wh/2 ~ '^h/i)- 
We found that the convergence rate in our case is very close to n = 2 for the scalar and Aq 
fields as expected for second order FDA, provided the numerical lattice is sufficiently dense to 
ensure we are in convergent regime. The convergence rate for for the Aq field was seen to be 
somewhat lower at, n ~ 0.7. Typical meshes that use to obtain results seen in|4]are of size 
Nw X Np = 400 X 800, which yields grid spacings of order hy, x hp 0.0025 x 0.1. 
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